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Abstract 

The modern theories of Grand Unification (GUT) and supersymmetric (SUSY) 
extensions of standard model (SM) suppose that the conservation laws of the SM 
may be violated to some small degree. The nuclei are well-suited as a laboratory to 
test fundamental symmetries and fundamental interactions like lepton flavor (LF) and 
lepton number (LN) conservation. A prominent role between experiments looking for 
LF and total LN violation play yet not observed processes of neutrinoless double beta 
decay (Oi^/3/3-decay). The GUT's and SUSY models offer a variety of mechanisms 
which allow Oz^/?/3-decay to occur. They are based on mixing of Majorana neutrinos 
and/or R-parity violation hypothesis. Although the Oi^/?/3-decay has not been seen it 
is possible to extract from the lower limits of the lifetime upper limits for the effective 
electron Majorana neutrino mass, effective right handed weak interaction parameters, 
the effective Majoron coupling constant, R-parity violating SUSY parameters etc. 

A condition for obtaining reliable limits for these fundamental quantities is that 
the nuclear matrix elements governing this process can be calculated correctly. The 
nuclear structure wave functions can be tested by calculating the two neutrino double 
beta decay (2z//3/3-decay) for which we have experimental data and not only lower 
limits as for the Oi^/3/5-decay. For open shell nuclei the method of choice has been the 
Quasiparticle Random Phase Approximation, which treats Fermion pairs as bosons. It 
has been found that by extending the QRPA including Fermion commutation relations 
better agreement with 2z//3/3-decay experiments is achieved. This increases also the 
reliability of conclusions from the upper limits on the Oi^/3/3-decay transition probability. 

In this work the limits on the LN violating parameters extracted from current 
decay experiments are listed. Studies in respect to future Oz//3/3-decay experimental 
projects are also presented. 

PACS numbers: 12.60.Jv, ll.30.Er, ll.30.Fs, 23.40.Bw 



1 Introduction 

The standard model (SM) represents the simplest and most economical theory which de- 
scribes jointly weak and electromagnetic interactions. It describes well all terrestrial exper- 
imental results known today. Nevertheless the SM cannot be considered as the ultimative 
theory of nature and is likely to describe the effective interaction at low energy of an underly- 
ing more fundamental theory. The supersymmetry (SUSY), which is one of the fundamental 
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new symmetries of nature, is believed to be next step beyond the successful SM. The super- 
symmetry is the symmetry between fermions and bosons, which has to be broken in order to 
explain the phenomenology of the elementary particles and their superpartners. It is the only 
known symmetry which can stabilize the elementary Higgs boson mass with respect to oth- 
erwise uncontrollable radiative corrections. The minimal supersymmetric model (MSSM), 
that leads to the SM at low energies has been the subject of extensive investigations. 

Both the SM and the MSSM leads to zero mass for the neutrinos. In view of the observed 
results on solar (Homestake Jll, Kamiokande [§, Gallex and SAGE 0]) and atmospheric 
neutrinos (1MB 0, Soudan 2 p, MACRO [0 and Super-Kamiokande 0) it is more appro- 
priate to consider the extensions of the SM and the MSSM that can lead to neutrino masses. 
Neutrino masses either require the existence of right-handed neutrinos or require violation of 
the lepton number (LN) so that Majorana masses are possible. So, one is forced to go beyond 
the minimal models again, whereby LF and/or LN violation can be allowed in the theory. 
A good candidate for such a theory is the left-right symmetric model of Grand Unification 
(GUT) inaugurated by Salam, Pati, Mohapatra and Senjanovic (especially models based 
on SO(IO) which have first been proposed by Fritzsch and Minkowski [|10[) and its super- 
symmetric version [jll|. The left-right symmetric models, representing generalization of the 



SU{2)l ® ^(1) SM, predict not only that the neutrino is a Majorana particle, that means it 
is up to a phase identical with its antiparticle, but automatically predict the neutrino has a 
mass and a weak right-handed interaction. The basic idea behind grand unified models is an 
extension of the local gauge invariance from quantum chromodynamics (SU3) involving only 
the colored quarks also to electrons and neutrinos. We note that the non-supersymmetric 
left-right models suffer from the hierarchy problem. 

The expectations arising from GUT's and theirs SUSY versions are that the conservation 
laws of the SM, e.g. LN conservation, may be violated to some small degree. In the left-right 
symmetric models the LN conservation is broken by the presence of the Majorana neutrino 
mass. The LN violation is also inbuilt in those SUSY theories where R-parity, defined as 
Rp = (_i)3B+-L+2S j^g^ g^j^j^ L are the spin, baryon and lepton number, respectively) is not 
a conserved quantity anymore. The conservation of LN is among the most stringently tested 
laws of physics nowadays. The nuclei are well-suited as laboratory to test this fundamental 
symmetry due to the fact that a variety of quantum numbers is available as initial and final 
states. The neutrinoless double beta decay (Oi//5/5-decay), 

{A,Z) ^ {A,Z + 2) + 2e-, (1) 

which involves the emission of two electrons and no neutrinos, has been long recognized as a 
powerfool tool to study the LN conservation. The Oz//3/3-decay takes place only if the neutrino 
is a Majorana particle with non-zero mass [|1^. The GUT's and R-parity violating SUSY 



models offer a plethora of the Oi//3/3-decay mechanisms triggered by exchange of neutrinos, 
neutralinos, gluinos, leptoquarks etc. [|13], If one assumes that one mechanism at a time 
dominates, the half-life of the Oz//3/9-decay can be written as 

iT^J,r' = \LNV\'J2P'''G''', (2) 

i 

where LNV is some effective LN violating parameter, P^'^ is the real part of the product of 
two nuclear matrix elements governing the Oz//3/3-decay and G^" is the integrated kinematical 
factor. The sum over i runs over different phase space integrals weighted by a corresponding 
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product of nuclear matrix elements. There are different LNV parameters, e.g. the effective 
electron-neutrino mass, effective right handed weak interaction parameters, effective Majoron 
coupling constant and R-parity violating SUSY parameters, which incorporate elements of 
the fundamental interaction of Majorana neutrinos and/or R-parity violating interaction of 
SUSY particles (see e.g. the recent review articles 0, Q). The value of these parameters 
can be determined in two ways: i) One can extract upper bound on the LNV parameter 
from the best presently available experimental lower limit on the half life of the 0i//3/3- 
decay T^^'^^^ after calculating the corresponding nuclear matrix elements, ii) One can 
use the phenomenological constraints imposed by other experiments, e.g. those looking for 
the neutrino oscillation, to evaluate the LNV parameter explicitly, which further can be 
compared with the extracted one. The Oi//3/3-decay constraints on LNV parameters must be 
taken into account by the theoreticians, when they build new theories of grand unification. 

At present the searches for Oz//5/5-decay are pursued actively for different nuclear isotopes, 
e.g. ^^Ge (Heidelberg-Moscow coll. p| and IGEX coll. |T6[), i°°Mo (NEMO 0], and 



ELEGANTS ^^^Cd (the INR exper. |2| and the NEMO exper. [0), ^^^Te and 

^^^Xe (the Gotthard Xe exper. [^). The sensitivity of a given isotope to the different LN 



violating signals is determined by the value of the corresponding nuclear matrix element 
connecting the ground state of the initial and final nuclei with J = O"*" and the value of the 
kinematical factor determined by the energy release for this process. In order to correctly 
interpret the results of Oz//3/3-decay experiment, i.e. to obtain qualitative answers for the LN 
violating parameters, the mechanism of nuclear transitions has to be understood. The Qv(3(3- 
decay nuclear systems of interest are medium and heavy open shell nuclei with complicated 
nuclear structure. To test our ability to evaluate the nuclear matrix elements that govern 
the decay rate, it is desirable to describe the two- neutrino double beta decay (2i//3/3-decay) 
allowed in the SM: 

(v4, Z) ^ (A, Z + 2) + 2e~ + 2z>. (3) 

We note that each mode of the double beta decay requires the construction of the same 
many-body nuclear structure wave functions. 

A variety of nuclear techniques have been used in attempts to calculate 2z//5/5-decay 
matrix elements, which have been reviewed recently in Refs. O, WM . Especially the Quasi- 



particle Random Phase Approximation (QRPA) and its extensions have been found to be 
powerful models, considering their simplicity, to describe nuclear matrix elements, which 
require the summation over a complete set of intermediate nuclear states. The recent 2vj3j3- 



decay calculations [^, |33, ^ including the schematic ones |139| , ^ manifest that the inclusion 
of the Pauli exclusion principle (PEP) in the QRPA improves the predictive power of the 
theory giving more reliable prediction of the 2vf3(3 decay probability. 

In this contribution we present the upper limits on some effective LN violating parameters 
extracted from the current experimental limits of the 0z//9/?-decay lifetime for A = 76, 82, 96, 
100, 116, 128, 130, 136 and 150 isotopes, which are quantities of fundamental importance. 
A discussion in respect to the sensitivity of a given Oz//5/3-decaying isotope to the different 
LN violating signals is presented. Some related nuclear physical aspects as well as studies 
in respect to future Oz//3/3-decay experiments are addressed. 



2 Neutrinoless double beta decay 
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Figure 1: Mechanisms for Oz//?/?-decay associated with the exchange of a Majorana neutrino: 
(a) the hght and the heavy neutrino mass mechanism, (b), (c) right handed current mecha- 
nisms, (d) the Majoron mechanism. The following notation is used: (^l(_r) and tL{K) 
are left- (right-) handed u-quark, d-quark and electron, respectively. W is vector boson (light 
or hypothetical heavy) and z/j (i=l,2...) is the Majorana neutrino. 



2.1 Majorana neutrino mixing mechanisms 

The presently favored models of grand unification are left-right symmetric models . They 
contain left- and hypothetical right-handed vector bosons Wf; and W^. The vector bosons 
mediating the left and right-handed interaction are mixed if the mass eigenstates and 
W2 are not identical with the weak eigenstates, which have a definite handedness.: 



cos C ■ Wf; sin ( 
- sin C ■ Wf; -t- cos C ■ 



w 



R 
± 
R ■ 



(4) 



C, is the mixing angle of the vector bosons. The left-right symmetry is broken since the vector 
bosons Wf" and W2 obtain different masses by the Higgs mechanism. Since we have not seen 
a right-handed weak interaction the mass of the heavy, mainly "right-handed" vector boson 
must be much larger than the mass of the light (81 GeV) vector boson, which is responsible 
for the left-handed force. 



4 



The weak interaction Hamiltonian must now be generalized. 



w 



Gp 

71 



{Jl ■ 3l) + tgajR ■ Jl) + tg^JL ■ Jr) + [jj^j {Jr ■ Jr) 



+ h.c, 



Jl/Jr = 1pplt,{gvT gAlb)i^n, ji/ji? = e7^(l =F 75) ^^eL/i?, (5) 

where gy = 1 and gA = 1.25. and are the hght and heavy vector boson masses, 

respectively. The capital Jl and Jr indicate the hadronic left- and right-handed currents 
changing a neutron into a proton, respectively. The lower case ji {jr) is the left (right) 
handed leptonic current which create an electron (or annihilate a positron) and annihilate 
left (right) handed current neutrino u^l {i-'er)- The weak interaction Hamiltonian (^ is 
given for C ^ 1 and ^ keeping only the lower order terms in expansion of tgC, and 
/M^ parameters. 

The left-right symmetric models allow us to explain the smallness of the neutrino mass 
within the so called see-saw mechanism in the most natural way. It is supposed that the 
neutrino mixing does take place according to 

^eL = E U^kXkL+ E U^kNtL, 
k=light k=heavy 

^eR = E U!},XkR+ E U!},N,n, (6) 

k=light k=heavy 

where, Xk {Nk) are fields of light (heavy) Majorana neutrinos with masses {rrik << 1 
MeV) and Mk {Mk » 100 GeV), respectively, and f/^^ are unitary mixing matrices. 
In the case of the most general lepton mixing originating from a Dirac-Majorana mass term 
in the Lagrangian the flavor neutrino fields are superposition of three light and three heavy 
Majorana neutrinos with definite mass. The fields Xk and Nk satisfy the Majorana condition: 
XkCk = C Xfc ? ^fc^fe = C , where C denotes the charge conjugation and ^, ^ are phase 
factors. 

The possible quark level neutrino mixing mechanisms of Oz//5/5-decay are displayed in 
Figs. ID (a), (b) and (c). If the Oz//3/3-decay is triggered by exchange of a light (heavy) 
left-handed Majorana neutrino [see Fig. || (a)] the corresponding amplitude of the process 
is proportional to the LN violating parameter < m^^ > {r]j^): 

light 

<m^> = E {^ek)^ rrik, 

k 

heavy 

V. = {uDHk^. (7) 

k k 

where rrip is the proton mass. The difference between < m^, > and rj^ comes from the 
fact that the neutrino propagator in the first and second case show different dependence on 



the mass of neutrinos |T^. We note that even if the neutrino is a Majorana particle but 
massless (i.e. there is no mixing of neutrinos) , the process in Fig. |1| (a) can not happen since 
for a pure left-handed weak interaction theory, the emitted neutrino must be right-handed 
(positive helicity), while the absorbed neutrino must be left-handed (negative helicity). With 
a finite mass the neutrino has not any more a good helicity and the interference term between 
the leading helicity and the small admixtures allows a Oi//3/3-decay. 
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The presence of the shght right-handed weak interaction in the GUT's allows the mech- 
anisms drawn in Figs. |1| (b) and (c). In this cases there is no helicity matching problem. 
The emitted neutrino from the left-handed vertex is right-handed as well as the absorbed 
neutrino at the right-handed vertex. Assuming only light neutrinos we distinguish two cases, 
i) The chiralities of the quark hadronic currents are the same as those of the leptonic currents 
coupled through the W-boson propagator [Fig. |l| (b)]. Thus the Oz//3/3-decay amplitude is 
proportional to 

/ A/tW \ 2 light 

<A> = (^^1 E U^Uf^i,. (8) 



k 



Recall that the W-boson propagator can be approximated by for Af = M|^,M^ [see 

Eq. ©]. 

ii) The chirality of the right-handed leptonic current is opposite to the coupled left-handed 
hadronic current [Fig. |1] (c)]. This is possible due to W-boson mixing. The corresponding 
effective LN violating parameter is 

light 

<V> = tgCY^ Ujl (9) 

k 

It is worthwhile to notice that the factor Y^lf^^ in Eqs. (||) and @ is expected 

to be small and even in the case there are only light neutrinos vanish due to orthogonality 



condition [Q . It indicates that the values of < A > and < rj > might be strongly suppressed 



assuming the see-saw neutrino mixing mechanism. 



2.2 Majoron mechanism 



The spontaneous breaking of the LN in the context of the see-saw model imply the existence 
of a physical Nambu-Goldstone boson [E3], called Majoron [HI], which is a light or massless 



boson with a very tiny coupling to neutrinos 

C^uu = T.^il5^j {i Irn 0) P,j, = ^ f/^ Ufp g^^p, (10) 

i<i a,f3=e,fi,T 

Here, Ui denotes both light Xi ^^id heavy A^, Majorana neutrinos. 

The Majoron might occur in the Majoron mode of the Oz//5/5-decay (Oi//?/?0-decay) 

{A, Z) ^ {A, Z + 2) + 2e- + (j). (11) 

and offers a new possibility for looking for a signal of new physics in double beta decay 
experiments. The Oj//3/30-decay mode yields a continuous electron spectrum for the sum 
of electron energies like the 2z//5/3-decay mode but differs from it by the position of the 
maximum as different numbers of light particles are present in the final state. We remind 
that in the case of Oi//3/3-decay a peak is expected to be at the end of the electron-electron 
coincidence spectra. 

The mechanism leading to a Oz//?/?(/)-decay mode is drawn in Fig. |l| (d). The experimental 
lower limits on the half-life of Oz//3/3(/)-decay allow to deduce the upper limit on the effective 
Majoron coupling constant < g >: 

light 

<9>=i: U^iU^jPij- (12) 



6 



2.3 R-parity violating SUSY mechanism 

Besides the simplest and the best known mechanism of LN violation based on the mixing 
of massive Majorana neutrinos advocated by different variants of the GUT's the R-parity 
violation proposed in the context of the MSSM is becoming the most popular scenario for 



LN violation (see e.g. reviews |i5, We remind that the R-parity symmetry assigns even 
R-parity to known particles of the SM and odd R-parity to their superpartners and that the 
Lagrangian of the MSSM conserves R-parity. The R-parity conservation is not required by 
gauge invariance or supersymmetry and might be broken at the Planck scale. The R-parity 
violation {Ijip) is introduced in the effective Lagrangian (or superpotential) of the MSSM in 
terms of a certain set of hidden sector fields. The trilinear part of the $,p superpotential 
takes the form 

= \,,UL,El + K^,UQ,Dl + KjiJJtDpi. (13) 

Here L and Q stand for lepton and quark doublet left-handed superfields while E'^,U'^, D'^ 
stand for lepton and up, down quark singlet superfields. Aj^fc, and X'l^^ are coupling 
constant and indices denote generations. The A"-terms are causing baryon number 

violation and the remaining ones the LN violation. In fact a combination of A' and A" leads 
to proton decay. 

If the SUSY models are correct the Oi//3/?-decay is feasible. The nuclear Oz//3/3-decay 
is triggered by the Oi//5/3-decay quark transition dd —>■ uu + e~e~. The relevant Feynman 
diagrams associated with gluino g and neutralino x trilinear SUSY contributions to the 
Oz//3/5-decay are drawn in Fig. ||. The ^p SUSY vertices are indicated with bold points. We 
see that the Ou (3 (3- decay amplitude is proportional to the X^^^ squared. 

There are two possibilities of the hadronization, i.e. coming from the quark level to 
the nucleon level. One can place the four quark into the two initial neutrons and two 
final protons, what is just the conventional two-nucleon mechanism of Oi//3/?-decay {nn 
pp + e~e~). This mechanism is strongly suppressed by the nucleon-nucleon repulsion at 
short distances for the exchange of heavy SUSY particles and heavy Majorana neutrinos. 
Another possibility is to incorporate quarks involved in pions in flight between nucleons. 



This possibility have been first pointed out by Pontecorvo [0. It is the so called pion 
exchange mechanism (vr" — 7r+ + e~e~). The pion-exchange mode leads to a long-range 
nuclear interaction, which is significantly less sensitive to short-hand correlations effects. It 
was found that ^pSUSY pion exchange contribution to the Oz//5/5-decay absolutely dominates 
over the conventional two nucleon mode realization P5| , pT| . 



The enhancement of the pion-exchange mode has also an origin in the bosonization of the 
7r~— > 7r+ -|- e~e~ vertex and is associated with the pseudoscalar JpJp and tensor Jj^ Jtiiv 
hadronic current structure of the effective R-parity violating Oi//3/3-decay Lagrangian on the 
quark level (Jp = ^756? and {Jj!" = ua'^''{l +'y^)d). The corresponding hadronic matrix 
elements are given as follows: 

5 

< 7r+(g) I JpJp 1 7r~(g) > ^ - < 7r+(g)|Jp|0 >< 0|Jp|7r-(g) > 

10 ,2 4 
9 {rriu + nidY 

< 7r+(g)|jrJTM.k"(g) > ~ -4<7r+(g)|Jp|0><0|Jp|7r-(g)>. (14) 
Here is the pion mass, /,r = 0.668 mj^. Taking the conventional values of the current 
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Figure 2: Feynman graphs for the supersymmetric contributions to Ou (3 (3-decay. ul, 
and cl have the same meaning as in Fig. [l|. ul, dji and are left-handed u-squark, 
right-handed d-squark and left-handed selectron, respectively, x 9 ai's neutralinos and 
gauginos, respectively. 



quark masses = 4.2 MeV, = 7.4 MeV one gets cp ~ 214. In the case of the exchange 
of a heavy Majorana neutrino there is a vector and axial-vector hadronic current structure 
JAV'^Av^l of the effective Lagrangian {J^y = u°''J^{i — l5)da)- We have 

< 7r+(g)|jyAMk-(g) > ^ ^<7r+(g)|Jp|0><0|Jp|7r-(g)> 

= = -KcAiq'). (15) 

Assuming the average momentum of the exchanged pion to be about 100 MeV we find 
ca ~ 0.61. We note that cp ^ ca- 

In order to derive a limit on the R-parity violating first generation Yukawa coupling X[ii 
from the observed absence of the Oz//3/3-decay it is necessary to use viable phenomenological 
assumptions about some of the fundamental parameters of the MSSM. In Ref. [|^) the 
ansatz of universal sparticle masses was assumed and that the lightest neutralino is bino-like. 
Within such phenomenological scenarios it was found that the gluino and neutralino exchange 



mechanisms are of comparable importance [13, 48]. Another possibility is to implement 
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relations among the weak scale values of all parameters entering the superpotential and the 
soft SUSY breaking Lagrangian and their values at the GUT scale. This scenario have been 
outlined in Refs. 



, 50 



It was shown that there is no unique answer to the problem of the 
dominance of neutralino and gluino contribution to Oz//?/3-decay. The dominance of any of 
them is bound with a different choice of the SUSY parameters mo and mi/2- It is worthwhile 
to notice that in the extraction of A'^j^^ the main uncertainty comes from the parameters of 



supersymmetry and not from the nuclear physics side |49, 51 



3 2z//3/3-decay and nuclear structure 

Since there are measurements available for the 2z//5/3-decay with the geochemical method 
{^^Se |31|], ^^Zr ^^^Te and ^^^Te H]) and with the radiochemical method {^^^U H 
and for seven nuclei even laboratory measurements ( '^^Ca [R3], '^^Ge pH|, ^"^Se [E9l, ^^Zr 



ri, 0, ig, ^o°Mo H, H, "^crf [H, g and ^^"^Nd H), one could try to calculate for 



a test of the theory the double beta-decay with two neutrinos and compare them with the 
data. We note that a positive evidence for a 2i//3/3-decay transition to the 0^ excited state 
of final nucleus was observed for ^°°Mo [^. 

The inverse half-life of the 2i//5/3-decay is free of unknown parameters on the particle 
physics side and can be expressed as a product of a phase-space factor G^'^ and the Gamow- 
Teller transition matrix element Mq^ in second order: 



[T^i%(Ot. - 0+ )]~^ = G'-^ \M^^\\ (16) 



where 



^ ^ <o;iA,(o)|i+><i+|A,(o)|o+> 

|0^ >, |0j > and |1^ > are respectively the wave functions of the initial, final and inter- 
mediate nuclei with corresponding energies Ei, Ef and En. A denotes the average energy 
A = {Ei — Ef)/2. Ak is the Gamow- Teller transition operator = Y.i^i~i^i)k, k=l,2,3. 

The calculation of Mgt remains to be challenging and attracts the specialists of different 
nuclear models. The computational complexity of Mgt consists in the reliable description 
of the complete set of the intermediate nuclear states. Recently, it has been shown the sum- 
mation over the intermediate nuclear states in the present 2z//3/5-decay studies corresponds 
to a summation over a class of meson exchange diagrams within the S-matrix approach |]52| . 

The nuclear shell model gives a satisfactory description only of the low lying excited 
states of nuclei. In the heavier nuclei there is a large number of basis states in the shell 
model which does not allow to perform realistic calculation without severe truncations. It 
is supposed to be the reason that the shell model predictions of 2z//3/3-decay rate for heavier 



nuclei, especially those for the Te region, show deviations from the experimental data |]53 |. 
We note that the feasibility of shell model calculations is growing with increasing computer 
facilities allowing to handle much larger configuration spaces. 

Many different nuclear structure aspects of the many-body Green function Mgt have 
been discussed. It was suggested by Abad et al. that Mgt could be dominated by 
the transitions through the lowest intermediate 1"^ state (so called Single-State-Dominance- 
Hypothesis (SSDH)). The SSDH could be realized in two ways:i) There is the true dominance 
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Figure 3: The upper part shows the way how in the Random Phase Approximation (RPA) 
the 2z//5/5 decay is calculated. For the Fermi transitions the (3~{n p) amplitude moves 
just a neutron into the same proton level and the (3~^{p —>■ n) amplitude moves a proton into 
the same neutron level. For the Gamow- Teller transitions it can also involve a spin flip, but 
the orbital part remains the same. One immediately realizes that the occupation and non- 
occupation amplitudes favor the (3~ amplitude, but disfavor the (3^ amplitude. There one 
has a transition from an unoccupied to an occupied single particle state, which is two-fold 
small (s^) first by the fact that the occupation amplitude for the proton Vp and secondly that 
the unoccupation amplitude for the neutron state are both small. Therefore the 2z//5/5 is 
drastically reduced. 



of the first 1^ state, i.e. the contribution from higher lying 1+ states to Mgt is negligible, ii) 
There is a cancelations among the contributions of higher lying l"*" states of the intermediate 
nucleus. The idea of SSDH have been outlined in Ref. showing that some experimental 



and theoretical evidence supports it for a few 2i//3/?-decay systems. 

The difficulty of the calculation of Mgt consists in the fact that the 2i//3/3-decay matrix 
elements are strongly suppressed. Its value is only small fraction of the double Gamow- Teller 
sum rule that scales roughly like the number of pairs of unpaired neutrons |]56| , 



Y,\< f\A- A\i>\'^ ^{N - Z){N - Z -I). (18) 
/ 

The proton- neutron QRPA (pn-QRPA) has been found successful in describing the suppres- 
sion mechanism for the 2z//3/3-decay |l57| , Fig. ^ explains why the 2i//3/3-decay amplitude 
is so drastically reduced. Therefore the small effects which normally do not play a major role 
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can affect the 2z//?/5-decay transition probability. If one looks to the second leg of the double 
beta-decay which is calculated backwards as a. l3~^ {p —>■ n) decay from the final nucleus to 
the intermediate nucleus one finds that the matrix elements involved in these diagrams are 
Pauli suppressed by a factor {unVpY = {smallY. The neutron-particle proton-hole force in 
the isovector channel, which is usually included is repulsive while the particle-particle force 
usually neglected is attractive. Therefore both excitations tend to cancel each other and 
therefore the amplitude [3^ is drastically reduced. This cancelation for the second leg could 
be even complete, i.e. the backgoing amphtudes and thus groundstate correlations cancel 
the leading forward going terms. 

Although one can obtain agreement within pn-QRPA with the measured 2vl3(3 data 
multiplying the particle-particle G-matrix elements of the nuclear Hamiltonian with a factor 
Qpp in a range of 0.8 < gpp < 1.2 ( gpp in principle should be equal to unity), two leaps of faith 
are usually quoted: i) The extreme sensitivity of Mgt to the strength of particle-particle 
interaction which does not allow a reliable prediction of the 2z//5/5-decay probability. We 
note that Mgt as a function of gpp crosses zero, ii) The collapse of the pn-QRPA solution 
within the physical range of gpp, what is supposed to be a phase transition. The collapse is 
caused by generation of too many ground state correlations with increasing strength of the 
attractive proton-neutron interaction. 

The study of the QRPA approximation scheme for different model spaces manifest that 
the problems i) and ii) are related with each other. The undesirable behavior of the 
pn-QRPA has its origin in the quasiboson approximation (QBA) violating the Pauli exclu- 
sion principle (PEP) and causing the QRPA excitation operators behave like bosons. The 
renormalized QRPA, which considers the PEP in an approximate way, shifts the collapse of 
the QRPA outside the physical range of gpp and shows a less sensitive dependence of M'qj, 
on gpp It allows us to predict more reliable values of the double beta decay matrix 

elements. The importance of the PEP for solving the problem of the QRPA collapse has 
been shown clearly within the schematic models, which are trying to simulate the realistic 



cases either by analytical solutions or by a minimal computational effort In Ref. 



iO| , to our knowledge for the first time, the solution of the QRPA equation with full consid- 
eration of the PEP was presented. It was found that restoring the PEP, the QRPA solutions 
are considerably stabilized and a better agreement with the exact solution is obtained. A 
new extension of the standard pn-QRPA "QRPA with PEP" was proposed, which consider 
the PEP in more appropriate way as the RQRPA and might work well also in the case of 
realistic calculations. 

We note that the calculation of the 2z//5/5-decay nuclear transition continues to be subject 
of interest, which stimulates the rapid development of the nuclear theory |^ |6^, 



4 Limits on LN violating parameters 

The limits deduced for LN violating parameters depend on the values of nuclear matrix 
element MEf'^, of the kinematical factor G^'^ and of the current experimental limit for a given 
isotope [see Eq. (H)]. Thus there is useful to introduce sensitivity parameters for a given 
isotope to the different LN violating parameters, which depend only on the characteristics 
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of a given nuclear system. There are the following: 



C<™.>(r) = 10^ |M^l^*>| v/Goi year, C,JY) = 10^ |M^XT| ^Goi year, 
C<x>{Y) = 10' IMO^^I v'Caa year, U,y{Y) = 10^ IM^^^I ^/C,, year, 
C<,>(r) = 10« |M^1^*>| v^Gb year, CAi^.l^-) = 10^ |M-^| v/^oi year. (19) 

The explicit form of M<^^*>, m'^^JL., ^gt-: C'aa, C*^?;, Cxx and Af^^ can be found e.g. in 
Refs. II, [1|, H. 

Admittedly there is a rather large spread between the calculated values of nuclear matrix 
elements within different nuclear theories (for ^^Ge the calculated rates differ by a factor 
of 3 [O, WW). In principle there are no exact criteria to decide which of them are correct. 
Nevertheless one can argue the RQRPA method offers more reliable results than the QRPA 
as the ground state correlation are better under control due to the consideration of PEP. 

The present limits on LN violating parameters < mj, >, rj^^, < X >, < rj > and < g > 
are associated with the two-nucleon mechanism for which the correct treatment of the weak 
nucleon current J'^^ is crucial. We have: 



9v{q^)l^ - igM{q^)^qu - gA{q^)l^l5 + gp{q^)q^l5 



^, (20) 



where = (p — p')^ is the momentum transferred from hadrons to leptons {p and p' are 
four momenta of neutron and proton, respectively) and cr^'-' = {i/2)[y^,'y'^]. gv{q'^), gAiiq"^), 
gA^q"^) and gp{q^) are the vector, weak-magnetism, axial-vector and induced pseudoscalar 
formfactors, respectively, which are real functions of a Lorentz scalar g^. 

The matrix elements M<m^*> and M<^!^> have been calculated by neglecting the role of 
induced nucleon currents (weak magnetism and induced pseudoscalar terms). Recently, it has 
been shown that they contribute significantly to the OuPfd-decay amplitude They modify 



the Gamow- Teller contribution and create a new tensor contribution. Their contribution is 
as important as that of unchanged Fermi matrix element. It was found that indeed such 
corrections cause a more or less uniform reduction of the M<^^^*> by approximately 30% 
throughout the periodic table. In the case of heavy neutrino exchange (M<^^>) the effect is 
much larger ( a factor of 3) [^]. The nucleon finite size has been taken into account through 
the phenomenological formfactors and the PCAC hypothesis. We note that in calculating the 
matrix elements involving the exchange of heavy neutrinos, the treatment of the short-range 
repulsion and nucleon finite size is crucial It is expected that the correct treatment 



of the induced pseudoscalar term (which is equivalent to a modification of the axial-vector 
current due to PCAC) might influence significantly also the amplitude of Oi//3/5-decay in the 
case of the right-handed current mechanisms. It goes without saying that the validity of this 
argument can be ultimatively assessed by numerical calculations. 

The numerical values of the sensitivity parameters Cx(^) =< >, Vn, < A >, 
< rj >, < g > and A'^^^^^) are listed in Table ^ In their calculation we used values of M<^j>, 
M<X>, and M^^ calculated within the pn-RQRPA |6|, ||. As there are no available 



pn-RQRPA results for Cxx, Cr^ri we used those of Ref. [^. The quantity Cx{Y) is an 
intrinsic characteristic of an isotope Y. The large numerical values of the sensitivity Cx{y) 
correspond to those isotopes within the group of decaying nuclei which are the most 
promising candidates for searching for the LN violating parameter X. However, we remind 
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that there are also other microscopic and macroscopic properties of the isotope, which are 
important for building a Oz//5/5-decay detector. By glancing the Table |I]we see that the most 
sensitive isotope is ^^'^Nd. It is mostly due to the large phase space integral and partial due 
to the larger nuclear matrix element We remark that the nucleus ^^^Nd is deformed 

and that in the calculation of the corresponding nuclear matrix element the effects of nuclear 
deformation, which might be important, were not taken into account. 

It is expected that the experimental constraints on the half-life of the Oz//3/3-decay are 
expected to be more stringent in future. Knowing the values of Cx(^) there is a straightfor- 
ward way to deduce limits on LN violating parameters X from the experimental half-lives 

rpOv—exp 
-^1/2 '■ 



< nil, > 
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10 



-5 
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:Ov—exp 
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<V>< 7 . -7^- 
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c<.>\ t; 
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pOu—exp 



(A' 



iiiJ 



1/2 
rriq 



10- 



10^4 years 



iOu~exp 



(21) 



K, is equal to 1.8 (gluino phenomenological scenario |]T3[). The normalization of 10 years 
was chosen so that the C's are of order unity. 

The current experimental upper bounds on the Oz//5/3-decay effective LN violating pa- 
rameters of interest for different isotopes are shown in Table |I]. We see that the Heidelberg- 
Moscow experiment |[T^ (In Table |l] we are giving the sensitivity of the experiment for ^^Ge 



being T^^ ^'^ ^ 1-6 x 10^^ as we want to compare this value with those from other exper- 
iments.) offers the most restrictive limit for < rrii, >, rj^, < X >, < r] > and the ^"^^Te 



experiment |^ for < g >. We note that if the experimental data from an other geochemical 
experiment on double beta decay of ^^^Te would be considered (T^/f = 1.5 x lO^^, see Refs. 
cited in [3), one would get less stringent limit on < g > {< g >< 1.5 x 10~^), which is 
comparable with the upper bound offered by the ^°°Mo experiment ||2T| (see Table [ip. 

At present the largest attention is paid to the effective electron Majorana neutrino mass 
parameter < m,y > in light of the positive signals in favor of oscillations of solar, atmospheric 
and terrestrial neutrinos. The masses and mixing angles can be determined from the available 
experimental data on neutrino oscillations and from astrophysical arguments by using some 
viable assumptions (hierarchical and non-hierarchical neutrino spectra etc ). At present the 
three and four neutrino mixing patterns are the most favorable ones [^, However, 
there is a discussion whether we really need the fourth sterile neutrino to fit the current 
experimental data. Knowing the elements of the neutrino mixing matrix one can draw 
conclusions about the Oi//3/5-decay, in particular on < m^, >, assuming neutrinos to be 
Majorana particles. From the study of different neutrino mixing schemes it follows that 
the upper bound on effective Majorana neutrino mass < m^, > could be ranging from 10"^ 
to 1 eV |7T|. From the Table |1| we see that the Heidelberg-Moscow Oz//3/3-decay 
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Table 1: The present state of the Majorana neutrino mass (hght and heavy), right-handed 
current, Majoron and I^p SUSY searches in Oz//?/3-decay experiments. Ty^~^^ and T^JP~'^'"f' 
are the best presently available lower limit on the half-life of the Oi>(3(3-decaj and 0i'j3f3(f)- 
decay for a given isotope, respectively. (x(X) denotes according to Eq. (19) the sensitivity 
of a given nucleus Y to the LN violating parameter X. The upper limits on < mj, >, 777V, 
<A>, <r7>, <(7> and A'^^^^^ are presented, gch.-geochemical data. 
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experiment []T5| implies < > to be less 0.5 eV. This fact allows us to make two important 
conclusions: i) The Oz//?/5-decay plays an important role in deciding among the alternative 
possibilities of neutrino mixing, ii) The lepton number violation is in reach of near future 
Oi//3/?-decay experiments, if the neutrino is a Majorana particle. An issue which only Oi//3/3- 
decay can decide. We remind that the discovery of the Oz//5/3-decay, what would be a major 
achievement for particle physics and cosmology, will implies only the upper bound on < rrii, > 
as a plethora of other Oz//5/5-decay mechanisms is in the game. It is supposed that only further 
measurements of Oz//3/5-decay transitions to the excited states of daughter nucleus together 
with inclusion of nuclear theory and study of different differential characteristics will allow 
us to decide which mechanism is the dominant one. 

The present neutrino experiments does not provide us with useful information about the 
mixing of heavy Majorana neutrinos (M^ >> 100 GeV). Therefore, it is difficult to extract 
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the mass of heavy neutrinos from the current hmit rj^ < 8.6 x 10~® (see Table |l|). If one 
assumes the corresponding Uei to be of the order of unity for the hghtest heavy neutrino 
mass Ml one gets: Mi > 1.1 x 10^ TeV. However, this element of the neutrino mixing 
matrix is expected to be rather small due to large differences in masses of light and heavy 
neutrinos within the see-saw mechanism. Therefore the real limit on Mk is supposed to be 
much weaker. It is worthwhile to notice that the limit on t]^ is extremely sensitive to the 
nucleon finite size and the short-range correlation effects |^ . The heavy neutrino exchange 



are 



nuclear matrix elements evaluated without inclusion of the induced nucleon currents ||T2 
considerably overestimated |Q and should not be used in deducing the limit on r]^. 

The present particle physics phenomenology does not allow us to deduce the magnitude 
of Y^k^^^ U^k f/J. entering the expressions for the effective right-handed current parameters 
< A > and < rj >. If we assume its value is about unity, then we get from the current limits 
on I < A > I < 7.5 X 10"'' and | < ?7 > | < 4.4 x 10~^ (see Table |T]) the corresponding bounds 
on the mass of the heavy vector boson and the mixing angle ( as follows: > 93. 
TeV and \tg(\ < 4.4 x 10~^. Mohapatra has shown that in the left-right symmetric model 
with spontaneous R-parity violation there is an upper limit on (in the limit tan C ^ 

M]y^ = MY) of at most 10 TeV By using this value for upper bound on (for the 



lower limit on we consider the value 100 GeV) we find 

light 



k 



U^k U^k 61 < 1.1 X 10-^ - 1.1 X 10-^ (22) 

\tgC\ < 3.8 X 10"^ - 3.8 x 10"^ (23) 

We remark that these limits could be modified after the pseudoscalar term of the nuclear 
current is properly taken into account. 

The Oz//3/3-decay offer the most stringent limit on the R-parity violating first generation 
Yukawa coupling A'^ |^. Its value depends on the masses of SUSY particles (see Eq. 
(|2ll) ). If the masses of squark m,q and gluino rrig would be at their present experimental 
lower bounds of 100 GeV we deduce from the observed absence of the Oz//3/5-decay A'^ < 
1.2 X 10""^ (phenomenological scenario). A conservative upper bound is obtained using the 
SUSY "naturalness" upper bound nig, ~ 1 TeV. It gives A'^ < 3.8 x 10~^ (see Table p. 
We mention that the limits on A'^ depend only weakly on the details of the nuclear structure 
as A'm is proportional to the inverse square root of the nuclear matrix element. In addition, 
the corresponding nuclear matrix elements are changing only slightly within the physical 
range of parameters of the nuclear Hamiltonian [HUl. However, A'^ depends quadratic on 
the masses of SUSY particles. In the GUT's constrained SUSY scenario there is a rather large 
SUSY parameter space. By using different sample of relevant SUSY parameters one ends 
up with significantly different limits on A'^ Finally, we stress that the above limits 

are very strong and lie beyond the reach of near future accelerator experiments (HERA, 
TEVATRON) [Q. However, we note that the collider experiments are potentially sensitive 
to other couplings A-^-^, X"jf^ etc. 

There are new experimental proposals for measurements of the Oi//3/3-decay for different 
isotopes. A new NEMO 3 experiment is under construction, which is expected to reach a 
lower limit on the Oz//5/3-decay half-life of the order of 10^^ years for ^"^Se, ^^Zr, ^°°Mo, ^^^Cd, 



^3°Te and ^^'^Nd nuclei in a period of about six years pi. The KAMLAND WB, CUORE 



76i and GENIUS experiments are under consideration at moment. The KAMLAND 



15 



Table 2: The expected limits on LN violating parameters of interest from the future 0^(3(3- 
decay experiments. The same notations as in Table I is used. 
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experiment is supposed to use as double beta decay emitter ^^^Xe isotope in a liquid scinti- 
lator and measure a half-life limit of about 5 x 10^^ years. In the CUORE experiment the 
cryogenic detector set-up will be made with crystals of Te02- The expected half-life limit 
for the Oz//?/3-decay of ^^°Te could reach the value 1 x 10^^ years. The largest half-life limit 
of 5.8 X 10^^ years may be achieved in the future GENIUS experiment by using one ton of 
enriched ^^Ge and one year for the measurement W^. If the above experiments could be 



realized one would get considerably stronger limits on Oz//3/3-decay lepton number violating 
parameters. They are listed in the Table ^ By glancing the Table ^ we see that NEMO 3 
^150^^-j and CUORE (^^°Te) experiments have a chance to reach the value of 0.1 eV {^^^Nd) 
for the effective neutrino mass and the GENIUS experiment could surpass this border to 
lower value of 0.03 eV. 



5 Summary 

In this contribution we have discussed the problem of lepton number violation in the context 
of rare nuclear processes, in particular of the Oi//3/?-decay, which has a broad potential for 
providing important information on modern particle physics. We have shown that the Ou^f]- 
decay has strong impact on physics beyond the Standard model in the way it constrains the 
parameters of other theories. The mechanism of LN violation within the Oi//3/5-decay has 
been discussed in the context of the problem of neutrino mixing and the R-parity violating 
SUSY extensions of the Standard model. The relevant LN violating parameters have been 
the effective Majorana neutrino mass < rrii, >, effective heavy neutrino mass parameter t]^, 
effective right-handed weak interaction parameters < A > and < t] >, effective Majoron 
coupling constant < g > and the first generation trilinear R-parity violating SUSY coupling 
Am. The restrictions on the lepton number violating parameters have be deduced from the 
current experimental constraints on Ou (3 (3-deca.Y half-life time for several isotopes of interest. 

The present limit on the effective Majorana neutrino mass < m,^ >< 0.5 eV deduced from 
the Oz//9/3-decay experiment have been discussed in connection with the different neutrino 
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mixing scenarios advocated by current data of neutrino experiments. We conclude that the 
0i//9/3-decay experiment plays an important role in the determination of the character of the 
neutrino mass spectrum. Some analysis in respect to the heavy Majorana neutrino have 
been also presented. By using the upper and the lower limit on the mass of the heavy vector 
boson constrained by the left-right symmetric model with spontaneous R-parity violation we 
have determined the allowed range for the mixing angle of vector bosons. It has been found 
that the current upper bound for the R-parity violating SUSY interaction constant X'm is 
< 1.2 X 10"'* (< 3.8 X 10~^) assuming the masses of SUSY particle to be on the scale of 100 
GeV (1 TeV). A discussion of the dominance of the pion-exchange R-parity violating mode 
for the Oz//3/3-decay process was also presented. 

The interpretation of the LN violating parameters involves some nuclear physics. It is 
necessary to explore the nuclear part of the Oi//3/3-decay probability. The predictive power 
of different nuclear wave functions can be tested in the 2z//5/5-decay. One needs a good 
description of the experimental data for the 2z//5/5 probability. We have discussed the recent 
progress in the field of the calculation of double beta decay matrix elements associated with 
the inclusion of the Pauli exclusion principle in the QRPA. The reliability of the calculated 
Oi//3/?-decay matrix elements was addressed. 

We found it useful to introduce sensitivity parameters CxiX) for a given isotope Y associ- 
ated with different LN violating signals, which are free of influences from particle physics and 
relate simply the experimental half-lives with LN violating parameters. The largest value of 
Cx{y) determines that isotope, which is the most sensitive to a given lepton number violating 
parameter X. It is an important information for planning new Oz//5/5-decay experiments. 

The Ot//3/3-decay offers with both theoretical and experimental investigations a view to 
physics beyond the SM. New Oz//3/5-decay experiments are in preparation or under considera- 
tion (NEMO 3, KAMLAND, CUORE, GENIUS). They could verify the validity of different 
mixing scheme of neutrinos. The expected limits on the LN violating parameters which could 
be reached in these experiments are presented in Table 0. However, there is a possibility that 
the Oz//3/3-decay could be detected in the forthcoming experiments. This would establish that 
the neutrino is a massive Majorana particle. The recent development in neutrino physics 
has triggered the hope that we could be close to this achievement. 
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